The Metal-Insulator Transition of NbC^: an Embedded Peierls Instability 
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Results of first principles augmented spherical wave electronic structure calculations for niobium 
dioxide are presented. Both metallic rutile and insulating low-temperature Nb02, which crystallizes 
in a distorted rutile structure, are correctly described within density functional theory and the local 
density approximation. Metallic conductivity is carried to equal amounts by metal t2 g orbitals, 
which fall into the one-dimensional du band and the isotropically dispersing bands. Hybridization 
of both types of bands is almost negligible outside narrow rods along the line X-R. In the low- 
temperature phase splitting of the d\\ band due to metal-metal dimerization as well as upshift of 
the bands due to increased p-d overlap remove the Fermi surface and open an optical band gap 
of about 0.1 eV. The metal-insulator transition arises as a Peierls instability of the dii band in an 
embedding background of electrons. This basic mechanism should also apply to VO2, where, 
however, electronic correlations are expected to play a greater role due to stronger localization of 
the 3d electrons. 

PACS numbers: 71.20.-b, 71.30.+h, 72.15.Nj 



Despite intense work over decades the metal-insulator 
transition (MIT) of VO2 has remained a matter of con- 
troversy ||, |, |, [| |, §, 0| . This is related to the simul- 
taneous occurrence of a structural transformation from 
the high-temperature rutile phase to a distorted mono- 
clinic structure, which is characterized by (i) pairing of 
the metal atoms within chains parallel to the rutile c axis 
and (ii) their lateral zigzag-like displacement S. Elec- 
tronic states near the Fermi energy are of mainly V 3d 
t-ig character. They separate into the d\\ band, which 
mediates V-V overlap along the metal chains, and the 
remaining bands At the transition, splitting of 
the dy band and upshift of the bands due to increased 
metal-oxygen overlap produce a finite band gap. Dis- 
pute centers about the question, whether the dn splitting 
is caused by metal dimerization or by increased electronic 
correlations resulting from the reduced screening by the 

electrons |3|, |ij . State of the art band calculations gave 
strong hints of a structural instability but were not able 
to reproduce the insulating gap due to the shortcomings 
of the local density approximation (LDA) S . 

Interestingly, only very few studies have dealt with the 
neighbouring dioxides, which display related phenomena. 
NbC>2 likewise undergoes a MIT (at 1081 K) and a simul- 
taneous structural transition from rutile to a distorted 
variant having a body-centered tetragonal (bet) lattice 
H [n|. Despite the differences in long range order, lo- 
cal deviations from rutile are the same as in VO2, i.e. 
niobium atoms dimerize and experience lateral displace- 
ments at the transition |yl [l^, In contrast, the 
metallic oxides M0O2, WO2, TcC>2, and a-Re02 all crys- 
tallize in the same monoclinic structure as VO2 PJ- No 
doubt the phase transitions of VO2 and NbC-2 as well 
as the destabilization of the rutile structure in all these 
dioxides call for a unified description ^ ||, [l0|. Yet, al- 
though work on the dioxides neighbouring VO2 would 
help clarifying the nature of the MIT, there are only few 



theoretical investigations in this direction |T^, |l6L p~7|| . 
For this reason, a complete and widely accepted picture 
of the rutile-related transition metal dioxides has not yet 
evolved. 

Here we present first principles electronic structure cal- 
culations for NbC>2 . These are the first calculations at all 
for the insulating phase and the first calculations for ru- 
tile NbC>2, which used the experimental crystal structure. 
Both phases are correctly described within density func- 
tional theory (DFT) and the LDA. The results strongly 
support the band theoretical point of view of the MIT as 
arising predominantly from a Peierls-type instability of 
the rutile phase. 

Evidence for the MIT being driven by strong metal- 
metal bonding was first inferred from the small and al- 
most temperature independent paramagnetic susceptibil- 
ity below the transition temperature as well as from the 
enhanced thermal expansion of the c lattice parameter 
just above T c §f, ||, |(], ||. Additional support came 
from strong decrease of T c , nonlinear variation of c, and 
crossover of the magnetic susceptibility to Curie- Weiss 
behaviour on substitution of Ti for Nb, which destroys 
the metal- metal bonding |p|, fo|, 

According to neutron-diffraction experiments the bet 
structure of insulating NbC>2 has 16 formula units per 
cell, space group I4i/a (CfiA and lattice constants 
a w 2\[2aji and c « 2cr [[ll], [ij, [lj|, where an = 
4.8464A and c R = 3.0316A are the rutile lattice param- 
eters |p2[ . The two different Nb-Nb distances of 2.7 and 
3.3A as arising from the metal-metal pairing along the 
rutile c axis |h| deviate considerably from the rutile 
value cr. In addition, the zigzag-like in-plane displace- 
ment of niobium atoms parallel to either the [110] or 
[110] direction of the rutile subcell leads to a variation 
of Nb-0 distances between 1.91 and 2.25A, which devi- 
ate from the high-temperature values of 2.00 and 2.08A 
p"2| , p2| . Critical scattering observed above T c at the 
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tetragonal wave vector qp = (1/4, 1/4, 1/2) pointed to a 
soft phonon mode, which, however, could not be clearly 
identified [jl3| [2^, |24j. Yet, using a shell model, Ger- 
vais and Kress were able to reveal softening of a P point 
phonon with a displacement pattern consistent with the 
low-temperature structure |25j . 

Room-temperature UPS and XPS experiments re- 
vealed a 9 eV wide occupied group of bands, which falls 
into the 1 eV wide Nb 4c? bands just below Ep and a 
6 eV wide group of O 2p bands at higher binding ener- 
gies |p6[ . The room temperature optical band gap was 
estimated to about 0.5 eV 0. By and large, photoe- 
mission data agree with band structure and Xa cluster 
calculations jlj], |l6[ ^7], |28), which, due to the complex- 
ity of the low-temperature structure, only addressed to 
the rutile phase. Only Sasaki et al. accounted for the 
bet structure by studying Nb 2 Oio clusters with either 
short or long Nb-Nb distance and deduced a band gap of 
0.68 eV; however, oxygen positions were not relaxed and 
the same value for all Nb-0 distances used |l6). Finally, 
calculated Fermi surfaces and generalized susceptibilities 
agreed with experiment in that they could not find a 
soft-mode instability Jl5], p8[ |. 

The present calculations were performed using the 
scalar-relativistic augmented spherical wave (ASW) 
method [^9|. Both the LDA and the generalized gradi- 
ent approximation were applied with essentially no differ- 
ences in the electronic structure. Crystallographic data 
as given by Bolzan et al. J2^] and Pynn et al. jl3| were 
used. In order to account for the loose packing of the 
crystal structures empty spheres, i.e. pseudo atoms with- 
out a nucleus, were included to model the correct shape 
of the crystal potential in large voids. Optimal empty 
sphere positions and radii of all spheres were automat- 
ically determined |50|. As a result, three and 14 in- 
equivalent empty sphere types with radii ranging from 
0.96 to 2.49 a_e were included in the rutile and bet cell, 
respectively, keeping the linear overlap of niobium and 
oxygen spheres below 18.5%. The basis set comprised 
Nb 5s, 5p, Ad and O 2s, 2p as well as empty sphere 
states. Brillouin zone sampling was done using an in- 
creased number of k points ranging from 18 to 1800 and 
12 to 552 points within the respective irreducible wedge. 

Partial densities of states (DOS) resulting from calcu- 
lations for the rutile and bet structures are displayed in 
Fig. ^ The gross features are similar in both phases. 
Three groups of bands are identified. While bands in 
the energy range from ss —9 to —3 eV are dominated by 
O 2p states, the two higher lying groups at and above 
the Fermi energy derive mainly from Nb Ad orbitals. All 
other states play only a negligible role in the energy inter- 
val shown. Crystal field splitting causes nearly complete 
energetical separation of the Nb id bands into two groups 
of t2 g and e y symmetry. Contributions of the O 2p and 
Nb Ad states to the upper two and the lower group, re- 
spectively, are indicative of covalent bonding. 
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FIG. 1: Partial DOS of (a) metallic and (b) insulating Nb0 2 
(per f.u.). Slight broadening is due to the DOS calculation 
scheme pi]]. 



Our results for insulating and metallic Nb02 are in 
good agreement with the photoemission experiments pff] 
and the previous calculations |l^, [H| [2?|, ^] . Differences 
with the latter could be traced back to the fact that these 
calculations were not self-consistent or, due to lack of 
crystallographic data, used a rutile structure obtained 
from symmetrizing the bet structure. 

The results for both phases differ with respect to the 
energetical position of the O 2p group of bands as well as 
the shape of the t^g partial DOS. In Fig. || these DOS are 
further separated into their symmetry components. To 
this end orbitals were transformed to a local coordinate 
system with the z and x axis parallel to the apical axis 
of the local octahedron and the rutile c axis, respectively 
fL7[ , While the d x 2_ y 2 (= d\\) and d yz states mediate 
er-type d-d overlap along the rutile c and a axes, respec- 
tively, the d xz states account for 7r-type overlap across 
the chains. Note that the latter two orbitals combine 
into the e y states. The t2 g groups of bands are displayed 
in Fig. ||, where the d x i_ y i bands are highlighted by bars 
with widths proportional to the d x i_ y i -contribution to 
the respective wave function. Thus, states with vanish- 
ing bars are of pure e y character (apart from small e y as 
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FIG. 2: Partial t2 g DOS of (a) metallic and (b) insulating 
NbO-2 (per f.u.). 



FIG. 3: Weighted band structure of (a) metallic and (b) in- 
sulating Nb02 (see text). 



well as O 2p contributions). 

Most striking in Fig. |](a) is the nearly perfect one- 
dimensional dispersion of the d x 2_ y i bands parallel to 
the rutile c axis, i.e. parallel to the niobium chains, which 
gives rise to the pronounced double peak structure in the 
partial DOS, Fig. ||(a). In contrast, dispersion of the e£ 
bands is of the same order in all directions. This leads to 
the single broad structure of the d xz partial DOS, while 
the d yz states fall into two maxima due to the aforemen- 
tioned inplane d-d overlap. The DOS at Ep is comprised 
from similar contributions from all three t2 g orbitals ex- 
plaining the rather isotropic electrical conductivity. Note 
that, except for the bands along the line X-R, the d\\ 
states hardly hybridize with the e g states and thus cou- 
pling between both types of bands is mainly via charge 
conservation. For this reason, the d x 2_ y 2 band may be re- 
garded as a one-dimensional band in a three-dimensional 
embedding background of e g bands. 

On going to insulating Nb02 two distinct changes oc- 
cur, i) The previously one-dimensional d x 2_ y 2 band is 
split into bonding and antibonding branches due to Nb- 
Nb dimcrization within the chains. As a consequence, 
the separation of the two peaks in the d x 2_ y i partial DOS 
has increased by more than 1 eV and the contributions 



in between have been considerably reduced, ii) The e g 
states experienced energetical upshift by « 0.5 eV due to 
increased p-d overlap, causing almost complete depopu- 
lation of these orbitals. This is due to the much reduced 
Nb-0 distances arising from the zigzag-like displacement 
of the Nb atoms. However, the lateral displacements also 
affect the inplane metal-metal bonding. For this reason, 
the two peak structure of the d yz partial DOS has van- 
ished and both e g partial DOS have become more similar. 
Taken together these changes lead to complete separation 
of the low lying d x 2_ y 2 bands from the e g states, which 
gives rise to a band gap of 0.1 eV. This value was de- 
termined from checking band extrema on a 16 x 16 x 16 
mesh within the Brillouin zone. Note that the hybridiza- 
tion between both types of bands is still negligible. 

To conclude, the electronic structures of both metal- 
lic and insulating Nb02 are well described within DFT 
and LDA. The near Ep electronic structure of metallic 
Nb02 consists of two very weakly hybridizing types of ti g 
bands, namely, d|| states with one-dimensional dispersion 
parallel to the Nb chains and isotropically dispersing e g 
bands. In the low-temperature structure Nb-Nb dimcr- 
ization splits the d|| band into bonding and antibonding 
branches, whereas the e g states shift to higher energies 
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due to reduced Nb-0 distances. Since coupling between 
both types of bands is still by charge conservation rather 
than hybridization, the insulating state is interpreted as 
due to a Peierls instability of the d» bands in an embed- 
ding reservoir of electrons. 

The results complement previous work on VO2 as well 
as M0O2 H and give rise to a unified picture for 
the early transition metal dioxides, which explains (i) 
destabilization of the rutile structure and (ii) the metal- 
insulator transition of the d 1 members. Since the 4d 
bands are broader than the 3d bands metal-metal bond- 
ing is stronger in the fourth row dioxides this leading to 
the higher transition temperature of NbC>2 as compared 
to VO2. In contrast, the absence of a gap in LDA cal- 
culations for the latter compound might be indicative of 
a considerable amount of electronic correlations. Here, a 
combination of LDA and dynamical mean field theory as 
recently successfully applied to V2O3 (32) might be more 
appropriate. 
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